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It was proposed recently by Murakami et al. [Science 301, 1348(2003)] that in a large class ol 
p-doped semiconductors, an applied electric field can drive a quantum dissipationless spin current 
in the direction perpendicular to the electric field. In this paper we investigate the eflects of spin 
imbalance on this intrinsic spin Hall effect. We show that in a real sample with boundaries, due 
to the presence of spin imbalance near the edges of the sample, the spin Hall conductivity is not 
a constant but a sensitively position-dependent quantity, and due to this fact, in order to take 
the effects of spin imbalance properly into account, a microscopic calculation of both the quantum 
dissipationless spin Hall current and the spin accumulation on an equal footing is thus required. 
Based on such a microscopic calculation, a detailed discussion of the effects of spin imbalance on 
the intrinsic spin Hall effect in thin slabs of p-doped semiconductors are presented. 
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I. INTRODUCTION 

Efficient injection and coherent control of spins in 
non-magnetic semiconductors represent two principal 
challenges in the emerging field of spintronics, a new 
paradigm of semiconductor electronics based on the uti- 
lization of the electron's spin degree of freedomi. At 
a first glance, it seems a trivial thing to inject spins 
into non-magnetic semiconductors by use of ferromag- 
netic metals as sources. However, in reality it is not 
practical because most of the spin polarizations will be 
lost at the interface between metal and semiconductor 
due to the large conductivity mismatch^'"^. A possible 
approach that can solve this problem is to replace ferro- 
magnetic metals by ferromagnetic semiconductors ( such 
as Gai-ajMn-cAs ) as sources of spin injectiorA^, but 
for practical use at room temperature, the Curie temper- 
atures of ferromagnetic semiconductors are still too low. 
Due to such difficulties, how to achieve efficient injection 
of spins into non-magnetic semiconductors at room tem- 
perature remains an open question and more great efforts 
are still needed. Recently, based on the Luttinger ef- 
fective Hamiltoniani, Murakami et al. theoretically pre- 
dicted that an extraordinary spin Hall effect may occur in 
a large class of p-doped semiconductors ( such as Si, Ge, 
and GaAs ), which means that in such a semiconductor, 
an applied electric field can drive a substantial amount of 
quantum dissipationless spin current in the direction per- 
pendicular to the electric field, and the spin current does 
not decrease substantially even at room temperature^.. 
This effect might reveal a new way for achieving efficient 
spin injection in non-magnetic semiconductors at room 
temperature without the need of ferromagnetic metals 
and may also find some other important applications in 
spintronics. Prior to the discovery of this effect, a sim- 
ilar effect was also predicted by Hirsch^" and discussed 
extensively by several other author From the the- 
oretical points of view, the effect conceived by Hirsch 



is an extrinsic spin Hall effect, which is caused by the 
spin-orbit dependent anisotropic scatterings from impu- 
rities but not an intrinsic property of a material, and 
it will disappear completely in the absence of impurity 
scatterings. The spin current generated by the extrinsic 
spin Hall effect was shown to be rather smal]pi2*ii*i^, so 
it is of little use in the problem of spin injections in non- 
magnetic semiconductors. Unlike the extrinsic spin Hall 
effect, the spin Hall effect proposed in Refs.[8-9] is purely 
intrinsic, which arises from the intrinsic spin-orbit cou- 
pling in the valence bands of p-doped semiconductors and 
does not rely on any spin-orbit dependent anisotropic 
scatterings from impurities. From a more profound point 
of view, this effect has a deep topological character and 
shares some basic features with the quantum Hall edge 
current both physically and mathematically^'^. For ex- 
ample, just like the case of quantum Hall cSeciM^i^^, 
the spin Hall conductivity due to this effect is a dissi- 
pationless transport coefficient and can be expressed as 
an integral over all states below the Fermi energy, and 
the contribution of each state can be expressed entirely 
in terms of the curvature of a gauge field in momentum 
spaceSiS. Due to such features, the spin current gener- 
ated by this intrinsic spin Hall effect can be very large ( 
comparable to the ordinary charge currents) and hence 
can serve as an effective source for efficient injections of 
spins in non-magnetic semiconductors at room tempera- 
ture. Very recently, a similar intrinsic spin Hall effect was 
also found by Sinova et al. in two-dimensional electron 
gases ( 2DEGs) with Rashba spin-orbit couplingi&. They 
found that in 2DEGs with Rashba spin-orbit coupling, 
the dissipationless and intrinsic spin Hall conductivity 
will take a universal value as long as both spin-orbit split 
bands are occupied. It is anticipated this effect will also 
find some important applications in the emerging field of 
spintronics. 

Although some basic concepts about the intrinsic spin 
Hall effect are cleari2i2ii&, there are still a number of 
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important questions which are needed to be further 
clarified, and in the last year many theoretical works 
have been devoted to the study of this extraordinary 
effect|i^ii8,i9i20,2ii22,23,24,25,26 Basically, most of these 

theoretical works have been focused on the calculation 
of the intrinsic spin Hall conductivity. In the present pa- 
per, we present a theoretical investigation on the effects 
of spin imbalance on the intrinsic spin Hall effect in p- 
doped semiconductors. While it was well known both 
experimentally and theoretically that in spin-polarized 
transport phcnomenEkSLSSiSSiSSiSi ( including in semicon- 
ductor spintronics devices^Si^^iM ) spin imbalance may 
have significant influences on the transports of spins, 
what influences spin imbalance will have on the intrin- 
sic spin Hall effect is still a new subject and has not yet 
been explored. For the intrinsic spin Hall effect, from 
both the experimental and theoretical points of view, 
a clear understanding of the effects of spin imbalance 
would be much desirable because spin imbalance may 
not only have some signiflcant influences on the electric- 
field driven quantum dissipationless spin current and on 
its practical applications but also play a crucial role in 
the experimental measurement of the effect^. In this pa- 
per, based on a solid microscopic ground, we will derive a 
set of self-consistent spin transport equations which will 
present a proper description on the interplay between the 
spin imbalance and the electric-field driven quantum dis- 
sipationless spin current in the intrinsic spin Hall effect in 
p-doped semiconductors. Starting from these spin trans- 
port equations and with the help of appropriate bound- 
ary conditions, the quantum dissipationless spin current 
and the induced nonequilibrium spin accumulation in an 
actual sample with boundaries can be calculated simul- 
taneously on an equal footing. Our results show that 
the characteristics of the interplay between the quantum 
dissipationless spin current and the spin imbalance in 
the intrinsic spin Hall effect in p-doped semiconductors 
are very different from what was found in usual spin- 
polarized transport phenomena ( including in the extrin- 
sic spin Hall effect ), and some usual concepts about the 
interplay between spin current and spin imbalance can- 
not be applied to the intrinsic spin Hall effect. 



The paper is organized as follows: In Sec. II, we will 
present a microscopic derivation of the spin transport 
equations for describing the intrinsic spin Hall effect 
in p-doped semiconductors. In our derivation, the ef- 
fects of spin imbalance will be included explicitly. In 
Sec. HI, by solving these spin transport equations with 
the help of appropriate boundary conditions, the electric- 
field driven quantum dissipationless spin current and the 
induced nonequilibrium spin accumulation in thin slabs 
of p-doped semiconductors will be calculated explicitly. 



II. SPIN TRANSPORT EQUATIONS IN THE 
PRESENCE OF SPIN IMBALANCE 



In a large class of p-doped semiconductors such as Si, 
Ge, and GaAs, the valence bands are fourfold degener- 
ate at the r point. In the momentum representation and 
taking the hole picture, the valance bands in such semi- 
conductors can be described by the following Luttinger 
effective HamiltonianL^ 



^o = ;^[(7i + x72)k'- 272(k-S)2], 
2m 2 
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where Si is the spin-3/2 matrix, and 72 are the Lut- 
tinger parameters. For a given wave vector k , the Hamil- 
tonian (1) has two eigenvalues, given by 



e/f(k) = eA=±3/2(k) 
eL(k) = eA=±i/2(k) 
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where X = h ■ S//c is a good quantum number of the 
Hamiltonian Hq- The hole bands described by Eqs.(2- 
3) are referred to as the light-hole (LH) and heavy-hole 
(HH) bands, respectively. When a uniform electric field 
E is applied, the full Hamiltonian will be given hy H — 
Hq + eE • X, where — e is the charge of an electron. The 
equation of motion for the light and heavy holes in a 
uniform electric field has been derived in much detail in 
Ref. [8], and in the semiclassical approximation ( i.e., the 
spin is treated as a classical variable and hence commutes 
with the current operator ), the following equation of 
motion was obtained therein, 
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where eiji is the usual fully antisymmetric tensor in three 
dimensions. The occurrence of the last term in Eq.(4) 
is unusual, it represents a "Lorentz force" in momentum 
space and is a natural generalization of the quantum Hall 
g£jfg(,^3ji4ji5 -^jij-gg dimensions^iS. It is just due to this 
"Lorentz force" in momentum space ( which makes the 
hole velocity noncoUinear with its momentum ) that the 
applied electric field will drive a quantum dissipation- 
less spin Hall current in the direction perpendicular to 
the electric field. From Eq.(4), one can get that in the 
low temperature limit and in the semiclassical approxi- 
mation, the net spin current due to both the LH and HH 
bands will be given by^ 
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where Jj denotes the net spin current following to the Xj 
direction with spin parallel to the Xi direction, ^^(k) is 
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the fiUing of holes in the band with helicity A, and cr^ is 
the spin Hall conductivity, which is given by 



127r- 



■(3fc^ 



(6) 



with kf^ and kf^ denoting the Fermi wave numbers in the 
LH and HH bands, respectively. In obtaining Eqs.(5-6), 
one has assumed that the fillings of holes in each band 
can be described by the simple Fermi-Dirac equilibrium 
distribution function. An alternative way of calculating 
the intrinsic spin Hall conductivity is by use of the Kubo 
formulci2iiSiiLi2iSii2S. Based on the Kubo formula, as was 
shown in Ref . [9] , the full quantum treatment of the non- 
commutativity between the quantum spin and current 
operator will lead to a quantum correction to Eq.(6). But 
if one takes the semiclassical limit, the result will become 
the same as was given by Eq.(6). 

Eqs.(5-6) are the central results of Refs.[8-9]. They are 
valid in the absence of spin imbalance. But in a real sam- 
ple with boundaries, when a spin current circulates in it, 
spin imbalance will be caused inevitably near the edge 
of the sample by the spin current, and in the presence 
of spin imbalance, the spin current may be significantly 
different from what was given by Eqs.(5-6), especially in 
the regions near the edges of the sample. The reason 
for this is that in the presence of spin imbalance, the 
fillings of holes in the LH and HH bands may deviate 
significantly from the corresponding cases in the equilib- 
rium state. In order to take the effects of spin imbalance 
properly into account, one should obtain the distribu- 
tion function strictly by solving the Boltzman transport 
equation, which describes the changes of the distribution 
function in a nonequilibrium state. In a nonequilibrium 
but steady state, the Boltzman equation reads 



X . V/a(x, k) + k . Vk/A(x, k) = - E(^)a- 
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where {d f \ / dt)^^'^^ is the collision term due to impurity 

scatterings, and x and k are the drift velocities of holes 
in the real space and in the momentum space, respec- 
tively. Similar to Ref. [8], in this paper we will confine 
our discussion to the semiclassical limit and weak exter- 
nal electric field ( i.e., in the linear response regime ) so 
that the semiclassical equation of motion given by Eq.(4) 

can be appliedSiS. The collision term {dfx/dt)^^^^J will 
be given by 



f dfx .{coll.) 



'^'^ «;«,(k,k>(6A(k)-6v(k')) 
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where w^^^^i (k, k ) is the probability of a hole to be scat- 
tered from the state |kA) into the state |k A ) due to 
impurity scatterings, and the impurity scatterings will 
be assumed to be isotropic and spin- independent. 



In the equilibrium state, the fillings of holes in each 
band are stable and can be described by the simple Fermi- 
Dirac equilibrium distribution function. When the exter- 
nal electric field is applied and the system turns into an 
nonequilibrium but steady state, the fillings of holes in 
each band will still be stable but different from what was 
described by the simple Fermi-Dirac equilibrium distri- 
bution function. In the presence of spin imbalance, the 
changes of the fillings of holes in the LH and HH bands 
will be caused primarily by two kinds of contributions. 
The first kind of contribution is caused by the drifts of 
holes in the external electric field, and the second kind of 
contribution is due to the occurrence of spin imbalance. 
In the linear response regime ( i.e., in a weak electric 
field ), the deviations of the fillings of holes in each band 
from the corresponding cases in the equilibrium state are 
small, and the two kinds of contributions will be inde- 
pendent and both be proportional to df^/dep ( here 
/o — {exp[/3(eA(k) — ep)] + 1}^^ is the usual Fermi- 
Dirac equilibrium distribution function with /3 denoting 
the inverse of temperature and e f the Fermi level in the 
equilibrium state ). Considering this fact and by use of 
the relaxation time approximation, in the linear response 
regime the nonequilibrium distribution function f\ (x, k) 
( in a nonequilibrium but steady state ) can be expressed 
as the following, 

a j-O a f 

/a(x, k) = /" + Ma(x)|^ + er[EA(x) • Va]|^, (9) 

where Va ~ Kk/mx is the velocity of holes; r is the 
total relaxation time of holes due to impurity-induced 
random scatterings; Ea(x) is the total effective field felt 
by a moving hole in the band with helicity A, which is the 
sum of the external electric field E and a band-dependent 
effective field induced by the spin imbalance in the sam- 
ple. The detailed definition of Ea(x) and t will be given 
below. The second term in Eq.(9) just characterizes the 
deviation of the filling of holes in the band with helicity A 
from the corresponding case in the equilibrium state due 
to the occurrence of spin imbalance in the sample, and 
the presence of this term is mathematically equivalent to 
introducing a band-dependent " shift " fix the Fermi 
level ep- (It should be noted that unlike the correspond- 
ing cases in usual spin-polarized transport phenomena, 
here fix does not relate directly to the spin accumulation 
because the label A does not correspond to a fixed spin- 
polarization direction in real space. ). The third term in 
Eq.(9) denotes the change of the filling due to the drifts 
of holes in the external electric field and in the presence 
of impurity scatterings. By inserting Eq.(9) into Eq.(7) 
and assuming that the impurity scatterings are isotropic, 
the Boltzman equation can be simplified to the following 
form, 

Va • [E + 1 V A^a(x) + t V (Ea(x) • Va)] 
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where , ' is a characteristic relaxation time defined by 
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which characterizes the probability ( given by r^^, ) for 
a hole in the band with helicity A to be scattered into 
the band with helicity A due to impurity scatterings. 
For simplicity, in the following we will assume that the 
intra-band-scattering relaxation time taa = ( inde- 
pendent of A ) and the inter-band-scattering relaxation 
time T^^'f^_^^-^ = T2 ( independent of A and A ), and as 
usual, the total scattering probability for a hole ( given 
by r^^, i.e., the inverse of the total relaxation time of a 
hole ) can be given by 



-1 -Vr-i, 
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which is assumed to be independent of the band label 
A. Multiplying both sides of Eq.(lO) by and then 
integrating both sides with respect to and with the 
help of Eq.(12), one can find that the total effective field 
felt by a moving hole in the band with helicity A should 
be given by 



Ea(x) =E+i va^aW 

e 



(13) 



Eq.(13) suggests that in the presence of spin imbalance, 
in addition to the external electric field E, conduction 
electrons will also feel an effective field proportional to 
the gradient of the band- and position-dependent shift 
in the Fermi level. After r and Ea(x) arc determined 
from Eqs. (12-13), the nonequilibrium distribution func- 
tion /a(x, k) will also be determined by Eq.(9). Then in 
the semiclassical limit the electric-field driven quantum 
dissipationless spin current can be obtained through the 
following formula, 
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where ij = hkj /m\+ejkiH'2X'^ -7 /2)kikk/k^ (see Eq.(4) 
) and Si^A(k) = h(Xki/k)/3 are the velocity and the spin 
of a hole with momentum k and helicity A, respectively. 
( Since the spin-3/2 matrix S in the Hamiltonian (1) is 
a summation of the spin angular momentum s with spin 
one-half and the atomic orbital angular momentum 1 with 
spin one, the expectation value of s should be one third of 
SM-). By substituting Eqs.(12-13) into Eq.(9) and then 
inserting Eq.(9) into Eq.(14), the following result can be 
obtained 



J-(x) = as{:>i)eijkEk, 



(15) 



where <Ts (x) is the spin Hall conductivity in the presence 
of spin imbalance, which is given by 



a«(x) = a" 



[3fciJ^^(x)-fcf/ii(x)]. (16) 



Here is the spin Hall conductivity in the absence of 
spin imbalance, which has been defined in Eq.(6), and 
A*//(x) = Ai3/2(x) + M-3/2(x) and ^^(x) = Mi/2(x) + 
A'-i/2(^)- Eqs. (15-16) show that the effects of spin im- 
balance on the quantum dissipationless spin current due 
to the intrinsic spin Hall effect in p-doped semiconduc- 
tors are very different from what was found in usual spin- 
polarized phenomena ( including the extrinsic spin Hall 
ggg(,^i0jiiji2 y First, in the presence of spin imbalance, 
the spin Hall conductivity due to the intrinsic spin Hall 
effect in p-doped semiconductors might not be a constant 
but a position-dependent quantity. This is a new feature 
that was not seen before. Second, for the intrinsic spin 
Hall effect in p-doped semiconductors, the change of the 
quantum dissipationless spin current due to the occur- 
rence of spin imbalance is determined directly by /i^(x) 
( i.e., the band-dependent " shifts " in the Fermi level ) 
but is independent of the gradients of /i^(x). This is also 
significantly different from what was found in usual spin- 
polarized transports ( including the extrinsic spin Hall 
effect). These unusual characteristics of the intrinsic spin 
Hall effect in p-doped semiconductors can be understood 
by the following arguments. According to Eq.(6), the 
spin Hall conductivity should be determined uniquely by 
the Fermi wave numbers kf^ and k^. In the presence of 
spin imbalance, because the spin imbalance will induce 
a position and band dependent shift fJ-xi^) in the Fermi 
level, the Fermi wave numbers kf^ and fc£ will also be 
position-dependent, and the changes of fcjj and kf^ due 
to the occurrence of spin imbalance will be determined 
directly by /i^(x). Due to this reason, in the presence 
of spin imbalance, the spin Hall conductivity will be a 
position-dependent quantity, and the change of the spin 
Hall current due to the occurrence of spin imbalance will 
be determined by (^) but independent of the gradients 
of iJ'xiy^)- Finally, it should be pointed out that because 
we have considered only isotropic and spinless impurity 
scattering, the mechanism of the generation of the spin 
Hall current described by Eqs. (15-16) is still purely in- 
trinsic, though there are some significant differences be- 
tween Eq.(16) and Eq.(6). In fact, one can check that in 
the linear response regime the impurity scattering term ( 
i.e., the third term ) in Eq.(9) does not contribute to the 
spin Hall conductivity given by Eq.(16). This point will 
be more clearly seen from the results presented in Sec. HI. 
Of course, if the impurity scatterings are spin-orbit de- 
pendent, then the total spin current will contain not only 
the intrinsic part but also contain an extrinsic part due 
to the spin-orbit dependent impurity scatterings through 
the mechanism proposed by Hirschi2iiiii^. 

In the ordinary charge Hall effect, the charge Hall cur- 
rent causes charge imbalance in a sample and results in 
charge accumulation. Similarly, in the spin Hall effect, 
the spin imbalance caused by the spin Hall current will 
result in nonequilibrium spin accumulation in a sample. 
Corresponding to the quantum spin Hall current given 
by Eqs. (15-16), the nonequilibrium spin accumulation in- 
duced by the quantum spin Hall current can be obtained 
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as the following, 



16ee|, dxj 



[CiMi(x)-3CHMH(x)],(17) 



e^T{k[f/67r^mL and C, 



H 



where Cl 

e^T(kfjY /QTT^m,H arc the ordinary charge conductivities 
of the light holes and the heavy holes, respectively. 
Eqs.(15-17) show that in the intrinsic spin Hall effect in 
p-doped semiconductors, both the quantum dissipation- 
less spin current and the spin accumulation will depend 
sensitively on /X;^(x), i.e., the band-dependent " shifts 
" in the Fermi level. To find out the equations that 
^ix{'x) should satisfy, one can substitute Eqs. (12-13) into 
Eq.(lO) and integrate both sides of Eq.(lO) with respect 
to k, then one will arrive at the following equation 

vVa W = ;^[4/Xa(x) - /za(x) - Mi(x)], (18) 



where D\ = vryr/S is a characteristic hole diffusion 
length and is the band-dependent Fermi velocity. In 
addition to Eq.(18), /X;^(x) should also satisfy the charge 
neutrality condition, which requires that the net changes 
of the charge density due to the flow of the quantum 
spin current, given by Sp^ = -e X);^ / ^5;^[/a(x, k) - 

f°{€x{k))], should be zero. This leads to the following 
equation 



Mh(x) = -(!!^)3/2^Jx). 



(19) 



Eqs.(15)-(19) are the central results of the present pa- 
per. They constitute a set of self-consistent equations 
from which both the quantum dissipationless spin cur- 
rent and the spin accumulation due to the intrinsic spin 
Hall effect in a real sample of p-doped semiconductors 
with boundaries can obtained simultaneously with the 
help of appropriate boundary conditions. 



III. INTRINSIC SPIN HALL EFFECT IN THIN 
SLABS OF p-DOPED SEMICONDUCTORS 

Eqs. (15-19) are rather general and in principle they 
can be applied to samples with any kind of geometries. 
In the experimental measurement of the Hall effect ( in- 
cluding the spin Hall effect ), a thin slab geometry ( i.e., 
the Hall bar ) is usually applied. In this section, starting 
from Eqs. (15-19), we will present a detailed theoretical 
investigation on the intrinsic spin Hall effect in a thin slab 
of p-dopcd semiconductors. Wc assume that the longi- 
tudinal direction of the slab is along the z axis and the 
transverse direction along the y axis and the normal of 
the surface along the x axis, respectively, and an external 
electric field is applied in the longitudinal direction 



of the slab. The thickness of the slab is assumed to be 
much smaller than the hole diffusion length D\ and the 
length of the slab is assumed to be much larger than the 
width, so that only spin current flowing to the y direc- 
tion ( i.e., in the transverse direction of the slab ) with 
spin parallel to the x direction need to be considered. 
The two boundaries of the slab are assumed to be lo- 
cated at y = ±w/2, and w is the width of the slab. In 
general, it is very difflciflt to solve Eqs. (15-19) analyti- 
cally. In order to get some explicit expressions for the 
spin Hall current and the spin accumulation, we assume 
that in Eqs. (15-19) the hole diffusion length D\ is A- 
independent ( i.e., Dx = D ) and niL — tuh- Then from 
Eqs.(18-19) one can see that Univ) = ^3/2(2/) + ^-3/2(2/) 
and ^^{y) = ^1/2(2/)) + 1/2(2/) can be expressed as 

t^Hiv) ^ -l^dy) = Aey'^'' + Be-y'^'', (20) 

where A and B are two constant coefficients that need to 

be determined by the appropriate boimdary condition. 
In this paper, we will consider the transverse open cir- 
cuit boundary condition. In the transverse open circuit 
boundary condition, the spin Hall current will b(^ zero at 
the two boundaries of the slab, i.e., Jy{y = ±'U7/2) = 0. 
Substituting Eq.(20) into Eqs.(15-16), the spin Hafl cur- 
rent Jy{y) can be expressed as a function of the coef- 
ficients A and B. Then by use of the transverse open 
circuit boundary condition, the coefficients A and B can 
be determined, and one can get that 



B 



2ef (3fc^ - fcf ) 



1 



3fc|J 



cosh{w/2D)' 



(21) 



After the coeflScients A and B are determined, the spin 
Hall current Jy{y) and the spin Hall conductivity as{y) 
will also be obtained by inserting Eq.(20) into Eqs. (15- 
16), and the results are given by 



Jyiy) = <^s{y)E. 
as{y) = a°[l - 



cosh{y/2D) 
cosh{w/4D)' 



(22) 
(23) 



Eqs. (22-23) show that, in the presence of spin imbalance, 
both the spin Hall current and the spin Hall conductivity 
might be highly position-dependent and might also de- 
pend sensitively on the hole diffusion length D and the 
width w of the sample. The spin Hall conductivity (Ts{y) 
will be maximum at the center of the sample ( i.e., at 
y = Q ) and tend to be zero at the edges of the sample. 
Two limiting cases will be especially interesting. The first 
case is that the hole diffusion length D is much larger 
than the width w of the sample. In this limiting case 
the spin Hall current will be very small, i.e., <Ts{y) — 
everywhere. The second interesting case is that w ^ D. 
In this limiting case, the maximum value of the spin Hall 
conductivity will be given by (Ts{y) — cr^ ( at y = ) 
and <Ts{y) as y ^ ±w/2. These features can be 
seen clearly from Fig.l, where wc have plotted the po- 
sition dependence of the spin Hall conductivity as {y) in 
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FIG. 1: Illustration of the position dependences of the spin 
Hall conductivity (Js{y) in three cases with different ratios of 
wjD. ( wjD = 50 for the solid line, wjD — 10 for the dashed 
line, and w/D = 1 for the dotted line. (Ja(y) is normalized 
by (7°, i.e., the spin Hall conductivity in the absence of spin 
imbalance.) 

three cases with different ratios of w/D. From Fig.l and 
Eq.(23), one can see clearly that if no boundaries exist ( 
i.e., w — > cx) and hence no spin imbalance occurs ), the 
spin Hall conductivity will be a constant and return to 
the same result as was given by Eq.(6), i.e., the spin Hall 
conductivity will not be m^Hg^i bylweak isotropic and 
spinless impurity scatteri-ngs. '^if? i^-in agreement with 
Ref. [8] and also in agreement with the result obtained by 
a more accurate calculation performed in Ref. [20]. It is 
interesting to note that recently a similar conclusion was 
also obtained for the intrinsic spin Hall effect in 2DEGs 
with Rashba spin-orbit coupling by both numerical sim- 
ulations and analytical calculations, which suggest that 
in the presence of weak ( isotropic and spinless ) impu- 
rity scatterings, the intrinsic spin Hall conductivity in a 
Rashba two-dimensional electron gas should still take a 
universal value, proving that the sample size exceeds the 
localization length^'^'^^. Of course, it should be pointed 
out that at present different views also exist on this prob- 
lem. For example, in Ref. [24] it was argued that the 
spin-orbit-coupling induced intrinsic spin Hall current in 
a Rashba two-dimensional electron gas should vanish in 
the presence of impurity scatterings, even if the impurity 
scatterings are weak and spinless. 

The quantum dissipationless spin current generated by 
the intrinsic spin Hall effect does not carry charges ( i.e., 
it is a pure spin current), so it is very difficult to mea- 
sure the quantum dissipationless spin current directly. 
An indirect but much more convenient way to detect the 
quantum dissipationless spin current is to measure the 



nonequilibrium spin accumulation induced by the quan- 
tum dissipationless spin current. The nonequilibrium 
spin accumulation induced by the quantum dissipation- 
less spin Hall current in a thin slab of p-doped semicon- 
ductors can be got by inserting Eqs. (20-21) into Eq.(17), 
and the following result can be obtained, 

in^rea^,E,{CL + iCn) s\nh{y/2D) 
"^^^ e^eF{kl + ?>kfj)Dcos\v{w/AD) ' ^ ' 

Eq. (24) shows that the spin accumulation will be linearly 
proportional to the spin Hall conductivity cr^ and also de- 
pend sensitively on the ordinary charge conductivities Cl 
and Ch of the light and heavy holes. It also have a sen- 
sitive dependence on the hole diffusion length D and the 
sample width w. According to Eq. (24), for a infinitely 
large sample without boundaries ( i.e., w — *■ oo ), no spin 
accumulation will appear ( i.e., Sx{y) — for any finite 
y ). This is different from what was found in a Rashba 
two-dimensional electron gas, where it was found that 
the application of an in-plane electric field would induce a 
homogeneous nonequilibrium spin accumulation without 
resort to the boundary effectsi^^*^Si2I From Eq.(24), one 
can see that for an actual sample with boundaries, the 
spatial distribution of the spin accumulation due to the 
intrinsic spin Hall effect would be highly inhomogeneous. 
The spin accumulation will be maximum at the edges of 
the sample and vanish near the center of the slab, and the 
spin accumulation at the edges of the sample will increase 
with the increase of the sample width w. This has been 
illustrated in Fig. 2. From Fig. 2 one can sec that if the 
sample width w is much smaller than the hole diffusion 
length D, the spin accumulation induced by the quan- 
tum spin Hall current will be very small. On the other 
hand, if the sample width w is much larger than the hole 
diffusion length Z?, the spin accumulation at the edges of 
the sample will be almost a constant, independent of the 
sample width. This will be a merit for the experimental 
measurement of the intrinsic spin Hall effect. To obtain a 
quantitative estimation on the order of the magnitude of 
the spin accumulation induced by the quantum dissipa- 
tionless spin Hall current in a real sample, let us consider 
some actual experimental parameters. We take the or- 
dinary conductivity Cl,h ^ lO^il^^cm^^ and the hole 
diffusion length D ~ lOnm and h/ep ~ Ifsec. These pa- 
rameters are typical of the holes in GaAs with the hole 
density n ^ IQ-'^^cm"'^. The width w of the sample is 
assumed to be 100//m ( much larger than the hole diffu- 
sion length ) and a current density jx ^ lO^A/cm^. By 
use of the parameters listed above, from Eq.(24) it can 
be estimated that the spin accumulation at the edges of 
the sample will be on the order of 10^^ — lO^^/i^cm"^. 
Such magnitudes should be large enough to be measured 
by some ordinary experimental methods, for example, by 
the method proposed in Refs.[8-ll]. Finally, it should be 
pointed out that a rough estimation of the spin accumula- 
tion due to the quantum dissipationless spin Hall current 
was also presented in the supporting online material for 
Ref. [8] based on a simple analysis by use of the usual spin 
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FIG. 2: Illustration of the changes of the spin accumu- 
lation at the edges of a sample with the variation of the 
sample width. ( The spin accumulation is normalized by 

diffusion equation, but there are some significant differ- 
ences between the results-^i^tained ip^hc present paper 
and the corresponding resfrts"^'e€Rrtc(/ therein. This can 
be seen by making a compari^S between Eq.(24) ob- 
tained in the present paper and Eq.(S16) presented in 
the supporting onhne material for Ref.[8]. For exam- 
ple, according to Eq.(24) obtained in the present paper, 
the spin accumulation will not only depend on the spin 
Hall conductivity but also depend sensitively on the or- 
dinary charge conductivities of the light and heavy holes; 
however, according to Eq.(S16) in the supporting online 
material for Ref.[8], the spin accumulation will only de- 
pend on the spin Hall conductivity but is independent of 
the ordinary charge conductivities of the light and heavy 
holes. Our results show that though the mechanism of 
the intrinsic spin Hall effect is purely intrinsic, impurity 
scatterings might have some significant influences on the 
effect in a real sample with boundaries, and the usual spin 
diffusion equation might not be very suitable for describ- 
ing this effect. In fact, from the microscopic calculation 
presented in Sec. II, one can see that in general the quan- 



tum dissipationless spin Hall current ( given by Eqs.(15- 
16) ) and the spin accumulation ( given by Eq.(17) ) due 
to the intrinsic spin Hall effect do not satisfy the usual 
spin diffusion equation. 

In conclusion, in this paper we have presented a de- 
tailed theoretical investigation on the effects of spin im- 
balance on the intrinsic spin Hall effect in p-doped semi- 
conductors. We have shown that in a real sample with 
boundaries, the spin Hall conductivity might not be a 
constant but a sensitively position-dependent quantity 
due to the occurrence of spin imbalance near the edges 
of the sample, and in order to take the effects of spin 
imbalance properly into account, a microscopic calcula- 
tion of both the quantum dissipationless spin current and 
the spin accumulation based on an equal footing is thus 
required. We stress that some usual concepts about the 
interplay between spin current and spin imbalance might 
not be suitable for describing the intrinsic spin Hall ef- 
fect. After some modifications, the theory presented in 
this paper might also be applied to investigate the ef- 
fects of spin imbalance in the intrinsic spin Hall effect 
in 2DEGs with Rashba spin-orbit coupling. Finally, it 
should be pointed out that though in the last year many 
theoretical works have been devoted to the study of the 
intrinsic spin Hall effect, many controversial issues still 
exist concerning some fundamental aspects of this ex- 
traordinary effect. Among them, a big controversial is- 
sue is that what is the correct definition of spin current in 
materials with intrinsic spin-orbit coupling^iiSiiSiSi. As 
was argued in Ref.[9] and in Ref.[19], there are some 
difficulties with the conventional definition of spin cur- 
rent in spin-orbit-coupled systems, but it seemed that 
up to now there are still no unanimous views about this 
questior&iSiiLiSiiSiSiiS^. ( In the present paper we have 
used the same definition of Refs.[8] ). Because no un- 
ambiguous experimental detections have ever been done, 
on the present stage such controversial issues are difficult 
to be clarified unambiguously by pure theoretical argu- 
ments. But it could be anticipated that by combining 
future experimental results with more accurate theoret- 
ical investigations, these controversial issues should be 
able to be clarified unambiguously in the near future. 
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